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The diffusion of a pulse of small grains in an horizontal rotating drum is studied through discrete 
elements methods simulations. We present a theoretical analysis of the diffusion process in a one- 
dimensional confined space in order to elucidate the effect of the confining end-plate of the drum. 
We then show that the diffusion is neither subdiffusive nor super diffusive but normal. This is 
demonstrated by rescaling the concentration profiles obtained at various stages and by studying the 
time evolution of the mean squared deviation. Finally we study the self- diffusion of both large and 
small grains and we show that it is normal and that the diffusion coefficient is independent of the 
grain size. 
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I. INTRODUCTION 
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One of the most surprising features of mixtures 
of grains of different size, shape or material is 
to segregate (i.e., to unmix) under 
of conditions 0, | i |. Axial 
been extensively studied 

M H UM and theoreti- 
I23,|23,l23. Yet, full understanding is 
still lacking. Axial segregation occurs in an horizontal 
rotating drum partially filled with an inhomogeneous 
mixture of grains. This phenomenon (also known as 
banding) is known to perturb industrial processes such 
as pebble grinding or powder mixing. After typically 
a hundred rotations of the drum the grains of a kind 
gather in well-defined regions along the axis of the drum, 
forming a regular pattern. When the medium consists 
of a binary mixture of two species of grains differing by 
their size, bands of small and large particles alternate 
along the axis of the drum. The bands of small grains 
can be connected through a radial core which runs 
throughout the whole drum J^, H?j . The formation of 
the radial core occurs quickly in the first few rotations 
as the small grains migrate under the surface. 

Unlike radial segregation, axial segregation requires a 
transport of grains along the axis of the drum. Since the 
grains are initially mixed, they must travel along the axis 
of the drum to form bands. This underlines the impor- 
tance of the transport mechanisms in a rotating drum. 
Much theoretical work as been devoted to axial segre- 
gation, most of which assumes a normal diffusion of the 
grains along the axis [H [H [H [H E . This hvpoth- 
esis has recently been challenged by Khan et al. 17] who 
have reported remarkable experimental results. These 
authors studied the diffusion of an initial pulse of small 
grains among larger grains in a long drum. Since direct 
visualisation is not possible (because the small grains are 
buried under the surface), these authors used a projection 
technique. Using translucent large grains and opaque 



small grains, they recorded the shadow obtained when a 
light source is placed behind the drum. They found the 
diffusion process to be subdiffusive and to scale approx- 
imatively as t 1 ! 3 . They also studied the self-diffusion of 
salt grains and found again a sub-diffusive process with 
an similar t 1 ! 3 power law. In this article, we report nu- 
merical findings on the diffusion of a pulse of small grains. 
Interestingly, our results are in contradiction with those 
of Khan et al. pjj since we observed normal diffusion. 



a) 



b) 




FIG. 1: (Color online) Snapshots of the DEM simulation 
taken after rotations (a), 30 rotations (b) and 100 rotations 
(c). A initial pulse of small (red oneline) grains is inserted a 
the center of a rotating drum otherwise partially filled with 
large grains (yellow oneline). For clarity we used here an ini- 
tial pulse of length I = 25ds in a drum of length L — 120ds- 
Note that longer drums are used later. 



The outline of the paper is as follow. First we will 
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describe the simulation method. A theoretical analysis 
of one-dimensional diffusion in a confined space is then 
presented. We present concentration profiles and mean 
squared deviation and show that the diffusion process 
is neither subdiffusive or super diffusive. An instability 
in the average position is described. Finally, we report 
results on the self-diffusion of both small and large grains. 



II. DESCRIPTION OF THE SIMULATION 

This article presents results based on the soft-sphere 
molecular dynamics (MD) method, one of the Discrete 
Elements Methods (DEM). This method deals with de- 
formable frictional spheres colliding with one another. 
Although not flawless, it has been widely used in the pas t 
two decades and has proven to be very reliable [20l l3(3| . 
Here we study the diffusion of a pulse of small grains 
among larger grains in an horizontal rotating drum (see 
Fig. 1). We follow the positions of individual small grains 
and in the following, it is understood that we constantly 
refer to the small grains. In particular, the "concentra- 
tion" means concentration of small grains. 



A. Parameters of our simulations 

The mixture consists of two species of ideally spherical 
grains differing by their size. The small grains have a 
diameter ds = 5 mm and the large grains have a diameter 
of 2ds- The density is the same for both kinds of beads: 
p = 0.6g/cm 3 . The length of the drum, L, is varied from 
L = 60ds to L = 420^5 and its radius is set to 20 ds- The 
rotation speed is set to 0.5 rots -1 . The grains are initially 
placed in a cubic grid. The small grains are placed in the 
middle of the drum around x = and the large grains fill 
up the space between the pulse of small grains and the 
end plates. The number of small grains can be varied but 
unless otherwise mentioned is set to N = 10 300, which 
corresponds to an initial pulse of length I = 25ds- The 
number of large grains depends of course on the length 
of the drum. The simulations run for typically a few 
hundreds of rotations. 

The rotation is started and the medium compacts to 
lead to an average filling fraction of 37%. Radial seg- 
regation occurs rapidly (typically after five rotations of 
the drum) and the initial pulse is buried under the sur- 
face. The concentration profiles can be obtained at any 
time step. Knowing the exact position of every grain al- 
lows one to accurately compute the average position (x) 
and the mean squared deviation (x 2 ). Note that peri- 
odic boundary conditions could be used in order to avoid 
the effect of the end plates. However, the confinement 
would still play a role since the position of individual 
grains would still be limited. Moreover it may lead to 
nonphysical spurious effects. 



B. The molecular dynamics method 

The forces schemes used are the dashpot-spring model 
for the normal force F™j and the regularized Coulomb 
solid friction law for the tangential force [3l| F\- : respec- 



tively, F% 



ij^j and F ij 
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where Sij is the virtual overlap between the two parti- 



cles in contact defined by: 5^ = Ri + Rj 



where 



Ri and Rj are the radii of the particles i and j and 
is the distance between them . The force acts when- 
ever is positive and its frictional component is ori- 
ented in the opposite direction of the sliding velocity. 
k™j is a spring constant, 7™- a viscosity coefficient pro- 
ducing inelasticity, \i a friction coefficient, 7* a regular- 
ization viscous parameter, and vfj is the sliding velocity 
of the contact. If kf- and 7™- are constant, the restitu- 
tion coefficient, e, depends on the species of the grains 
colliding. In order to keep e constant the values of fcg- 
and 7™- are normalized using the effective radius R q q 
defined by 1/R e g = l/Ri + 1/Rji fcg = k%R /R e Q and 
lij ~ 7o ^gff/^o- The particle/ wall collisions are treated 
in the same fashion as particle/particle collisions, but 
with one particle having infinite mass and radius. The 
following values are used: Ro = 4 mm, kfi = 400 Nm" 1 , 
(leading to e - 0.9), 7* = 6 kgs -1 and \i = 0.3. The 
value of e was varied (from 0.4 to 0.9), which seemed to 
have only very little influence on the diffusion process. 

The equations of motion are integrated using the Verlet 
method with a time step dt = 1/30 At, where At is the 
duration of a collision (At « 10 _3 s). The simulations are 
typically run for 10 7 time steps, corresponding to a few 
hundreds of rotations. 



III. ID CONFINED DIFFUSION 

Before analyzing any results, one should quantify the 
influence of the end-plates. The limited space imposes 
some constraints on the diffusion process: the position x 
along the axis of the drum can only range from —L/2 to 
L/2 which imposes a limit to the mean squared deviation. 
In order to clarify the effect of the confinement we present 
in this section a theoretical analysis of an initial Dirac 
distribution (corresponding to the initial pulse) diffusing 
in a confined ID space. Note that an initial pulse function 
could also be used. However, this would add a degree 
of complexity to the problem without much benefit. In 
this brief section we only intend to qualitatively describe 
the effect of the confinement rather than studying it in 
details. The concentration, defined for x G [—L/2, L/2], 
is given by [32^ : 



^ 00 



-(x + nL) 2 
exp 4L>t 



n= — 00 



Note that the exact solution for an initial pulse of nonzero 
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width can be obtained by convoluating this solution to 
the initial pulse function. 




FIG. 2: (Color online) Left: concentration profiles obtained 
at various stages for a confined ID diffusion process. Right: 
mean squared deviation versus dimensionless time. At short 
times the function is a linear function of time but at longer 
time (x 2 ) saturates. Inset: the data plotted on a log- log scale. 




x/L 

FIG. 3: (Color online) Concentration profiles at different time 
obtained with a drum of length L = 420 ds . 



We can now plot the concentration profiles for various 
times (see Fig.|2j). Note that the time and position can be 
made dimensionless (using t* = tL/D 2 and = x/L). 
One can see that at short time the distribution is almost 
a Gaussian but at longer times, the confinement plays 
a important role. The distribution flattens and eventu- 
ally reaches a constant value of c(x) = 1/L, which corre- 
sponds to a mean squared deviation (x 2 )^ = L 2 /12. Fig- 
ure El also shows the time evolution of the mean squared 
deviation defined here by (x 2 ) = J Q c(x,t)x 2 dx. The 
function is clearly linear at short time (as confirmed by 
the log-log scale inset) but saturates a longer times. This 
behavior was expected but can lead to erroneous conclu- 
sions. Indeed, the curvature observed in (x 2 )(t) could be 
mistaken for a sign of a subdiffusive process whereas it is 
solely due to the confinement. Note that Khan et al [T^J 
used a rather long tube in their experiments, meaning 
that the confinement did not affect their results. 



IV. RESULTS 
A. Concentration profiles 

We now present results obtained in our 3D numerical 
simulations. The position of every grain at any time step 
is known, which allows one to compute the concentration 
profiles. The drum is divided in virtual vertical slices of 
length ds and the number of grains whose center is in 
a given slice is computed. Figure 03 shows concentration 
profiles of small grains measured at short times for a 
drum of length L = 420 ds- Note that the profiles tend 
toward Gaussian distributions as the grains slowly diffuse 
in the drum. 

Note that no axial segregation is visible at any point. 
Indeed, except for the obvious one at x = 0, there is no 
peak in the concentration. We believe that the number 
of small grains is too small compared to that of large 
grains to trigger axial segregation. 



A very good test to check whether the diffusion pro- 
cess is subdiffusive, normal or superdiffusive is to try to 
collapse various concentration profiles at different times 
onto one unique curve. This should be done by divid- 
ing the position by some power of the time (t a ) and for 
reasons of normalization (i.e. mass conservation) by si- 
multaneously multiplying the concentration by t a . The 
process is subdiffusive if a < 1/2, normal if a = 1/2 and 
super-diffusive if a > 1/2. Using the profiles taken at 
"short times" (i.e., before the confinement plays a role), 
we were able to obtain an excellent collapse of our data 
using the value a — 1/2. Figures Ufa) and ^b) shows 
the data in a linear and semilog scales. Figures Etc) and 
QJd) shows similar plots obtained using a = 1/3, as sug- 




FIG. 4: (Color online) Rescaled concentration profiles (in 
arbitrary units) at different times obtained with a drum of 
length L = 420 ds- The data is rescaled using t 1//2 for (a) 
and (b), and t 1 ^ 3 for (c) and (d). (a) and (c): linear scale, 
(b) and (d): semilog scale. The solid line on (a) and (b) is a 
Gaussian fit at t — 225 rotations. On can see that the t 1//2 
rescaling leads to a good collapse unlike the t 1 ^ 3 rescaling. 
This is particularly clear on the semilog plots. 
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gested by Khan et al. 17]. The first conclusion that can 
be drawn is that the collapse is excellent for a = 1/2 
and poor for a = 1/3. This is particularly clear on the 
semilog plots. The profiles on Fig. Sfd) tend to nar- 
row with increasing time whereas no trend is visible on 
Fig.Efb). This clearly shows that in our simulations the 
diffusion is a normal diffusion process. Moreover, the col- 
lapsed data can be well fitted by a Gaussian distribution, 
as clearly evidenced by the semilog plot on Fig.^b). 



B. Mean squared deviation 

The collapse of the different curves on Fig. B] is a good 
indication that the diffusion process in normal. An even 
better test is to plot the time evolution of the mean 
squared deviation defined here by (x 2 ) = 1/N — 
x®) 2 where is the initial position of the grain num- 
ber i. Figure [SK is a plot of the mean squared deviation 
versus time obtained for various values of L. Of course 
the saturation value (x 2 )^ is different for different val- 
ues of L but one can see that the initial slope is the same 
for all curves. This shows that the diffusion coefficient 
is well defined at "short times" and does not depend on 
the length of the drum. Figure EJb) is a rescaled plot of 
the same data. The time is divided by L 2 and the mean 
squared deviation by (x 2 )^. One can see that all data 
collapse at short times, showing again that the diffusion 
coefficient is well defined and that it is independent of L. 

One can see on Fig. |5Jb) that the mean squared 
deviation does not reach the theoretical saturation value 
(x 2 )oo. This is an indication that in the final state 
the concentration is not uniform. Indeed although no 
segregation bands are visible there is still a signature 
of axial segregation. The concentration profile in the 
steady state obtained for L = 60 ds is shown on Fig. 
One can see that the initial central peak has diffused 
but the concentration is not even in the final state: the 
concentration drops near the end plate, which has also 
been observed experimentally. This explains why the 
saturation value (x 2 )^ = L 2 /12 is not reached in our 
simulations. Note, however, that for longer drums this 
effect becomes negligible. 

The rescaled plot on Fig.JSJ} allows one to decide which 
time frame should be considered "short times". In par- 
ticular it is obvious from that picture that the run per- 
formed with L = 420 ds has not yet been affected by 
the end-plates after 200 rotations. Therefore, it allows 
one to study the diffusion process without considering 
the consequences of the confinement. Figure Q shows the 
mean squared deviation versus time in linear and log-log 
scales for a drum length L = 420 ds- Both plots clearly 
demonstrate that the diffusion is normal. Indeed, the 
data on Fig. [7] is well- fitted by a straight line on a linear 
scale. Similarly, after a short transient corresponding to 
the time it takes for radial segregation to be completed 
(approximatively 10 rotations here), the data plotted in 
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FIG. 5: (Color online) (a) Mean squared deviation versus 
time (in arbitrary units), (b) Rescaled data. The different 
curves correspond to different values of the drum length: Blue 
oneline: L = 60 ds, green online L = 150 ds red oneline : 
L — 420 ds- The solid line in (b) is the theoretical prediction. 




FIG. 6: (Color online) Steady state concentration profile for 
L = 60ds after 150 rotations. The distribution is not perfectly 
uniform: the concentration in small grains drops near the end 
plates. 



a log-log scale is perfectly fitted by a line of slope 1 for 
over a decade. These observations show without a doubt 
that the diffusion process in our simulations is neither 
subdiffusive nor superdiffusive but simply normal. 



5 




t (rotations) t (rotations) 

FIG. 7: (Color online) Mean squared deviation versus time in 
a linear (right) and log- log (left) scale for L — 420ds- Except 
for a brief transient during which radial segregation occurs 
(x 2 ) is a linear function of time the data plotted in a log- log 
plot scale is extremely well fitted by a line of slope 1. 



V. SELF-DIFFUSION 

In this section we will study the self-diffusion process. 
We would like to compare results obtained with three 
runs. The first run consists of monodisperse small 
particles, the second of large particles and the third one 
is the pulse experiment described above. The rotation 
speed, filling ratio and drum size are identical all three 
cases. In order to save computational time, we used a 
rather short drum (L — 60ds). The major difference 
between a self-diffusion numerical experiment and the 
diffusion of a pulse is that no radial segregation can exist 
in the self-diffusion. Since the mixture is monodisperse 
there cannot exist any form of segregation. It is therefore 
interesting to compare the two numerical experiments 
and elucidate the role of the radial segregation in the 
diffusion process. 

What should be measured in a self-diffusion experi- 
ment? Of course, there is no initial pulse in the drum but 
one can be arbitrarily defined. Indeed, one can track the 
grains whose initial position were within a given distance 
from the center of the drum. This would allow to plot 
"virtual concentration profiles." More simply, we can use 
the definition used before for the mean squared deviation 
(x 2 ) = 1/N — x®) 2 . Note, however, that the end 

plates play the same role as before. Indeed, the position 
is still limited. Moreover, the grains initially located near 
an end plate will "feel" the end plates at early stages. It 
is therefore necessary to measure the mean squared de- 
viation of grains initially centered around x = 0. For 
both self-diffusion simulations we used the grains located 
within a distance of 10ds from the center of the drum 
(x = 0). 

Figure [SI is a plot of {x 2 ) versus time for all three runs. 
Rather surprisingly, all three curve collapse. This shows 
that the diffusion process is independent of the grain size, 
which is very interesting. One could expect the diffusion 
coefficient to scale with the grain diameter and the two 



self-diffusion coefficients to be different but it is obviously 
not the case. Instead, our simulations show that ds is not 
a relevant length scale regarding the diffusion process. 

Maybe even more interestingly, the data collapse shows 
that radial segregation does not seem to affect the diffu- 
sion process. Indeed, the diffusivity is the same in the 
pulse experiment (where radial segregation exist) as in 
the self-diffusion experiments (where no segregation can 
exist). 
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FIG. 8: (Color online) Mean squared deviation versus time 
for the pulse experiment (blue oneline), self- diffusion of small 
(red oneline) and large (black) grains. One can see all three 
curves collapse, showing that the self- diffusion is normal and 
that the diffusion coefficient is independent of the grain size. 



VI. DISCUSSION AND CONCLUSION 

The diffusion of an initial pulse seems to be very helpful 
in understanding the basic mechanisms leading to axial 
segregation. It is however difficult to obtain precise mea- 
surements in an experimental system. Khan et al. [l7^ 
developed a clever projection technique which allows one 
to observe the hidden core of small grains. The shadow 
projected by the opaque small grains is clearly related to 
the concentration in small grains but one can question 
whether the link between the two is a linear relation or 
a more complex one. 

Our numerical simulations show consistent results in- 
dicating that the diffusion is normal. This conclusion is 
strongly supported by a collapse of concentrations pro- 
files rescaled by t 1 ! 2 . Moreover, the mean squared devi- 
ation is clearly a linear function of time as long as the 
confinement does not play a role. These results are in 
strong contradiction with those obtained experimentally 
by Khan et al. [l7|- These authors found a subdiffusive 
process which scales approximatively as t 1 / 3 . The origin 
of this discrepancy remains unknown. 

It could originate from the force models used in our 
simulations. In particular, it would be very interesting 
to test the consistency of our results using the Hertz mod- 
elor a tangential spring friction model |30( . The discrep- 
ancy might also originate in the projection method used 
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FIG. 9: View of the small grains projected perpendicularly to 
the free surface for L — 280(4 after 75 rotations (i.e., before 
the confinement plays a role). 



by Khan et al. [Tjj which is an indirect measurement of 
the concentration. However, it seems unlikely that the 
details of the projection techniuqe would change the dif- 
fusion power laws. One interesting check for that method 
would be to apply the projection technique to our numer- 
ical data. Figure EH shows a projection (perpendicularly 
to the free surface) of the core of small grains obtained 
for L — 280<i s after 75 rotations. This image looks rather 
differents from those of Khan et al. ^tJ- I n particular, 
one can observe individual grains, which is not possible 
using their projection technique. Applying the projec- 
tion technique to our numerical data is not straightfor- 
ward and the definition of such a procedure would obvi- 
ously drastically influence the outcome. Note, however, 
that Khan et al. have also observed a subdifTusive self- 
diffusion using direct imaging, which does not involve 
projection. This supports the idea that the 1/3 power 
law which they reported is not an artifact of their pro- 
jection technique and that the origin of the discrepancy 
between the experiments and the simulations may be due 
to some unaccounted for physical effect. 

After stimulating discussions with these authors, we 
believe that there are a number a physical differences 



between our two system that might lead to different 
results. One difference between the two systems is 
the ratio of the drum diameter to the avarage particle 
diameter S = D/d. In our simulations S — 26 and in 
their experiments S « 100. The discrepancy might also 
be caused by the difference in the shape of the grains 
although a subdifussive process was observed experi- 
mentally using bronze beads. Finally, the experiments 
were conducted in a humidity-controlled room and the 
capillary bridges between the grains might have modified 
the diffusion process. 

In conclusion, our numerical results showed that the 
diffusion along the axis of the drum is normal. Having 
studied the effect of the confinement, we can conclude 
that the curvature observed in the mean squared devi- 
ation is not a sign of subdiffusion. The rescaled con- 
centration profiles lead to the same conclusion that the 
diffusion process is normal. The study of self-diffusion 
shows that the diffusivity is independent of the grain size 
and is not affected by radial segregation. This results 
may shed some light on the mechanisms of formation of 
segregation bands since it indicates that the transport of 
grains along the axis of the drum is identical for both 
species of grains. We hope that these results will inspire 
theoretical models and help understanding the puzzling 
phenomenon of axial segregation. 
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